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Existence of weak solutions up to collision 
for viscous fluid-solid systems with slip 

David Gerard- Varet* Matthieu Hillairet^ 
July 3, 2012 

Abstract 

We study in this paper the movement of a rigid soUd inside an incompressible Navier-Stokes 
flow, within a bounded domain. We consider the case where sHp is allowed at the fluid/solid 
interface, through a Navier condition. Taking into account slip at the interface is very natural 
within this model, as classical no-slip conditions lead to unrealistic collisional behavior between 
the solid and the domain boundary. We prove for this model existence of weak solutions of 
Leray type, up to collision, in three dimensions. The key point is that, due to the shp condition, 
(~| . the velocity field is discontinuous across the fluid/solid interface. This prevents from obtaining 

global bounds on the velocity, which makes many aspects of the theory of weak solutions for 
Dirichlet conditions unadapted. 

1 Introduction 

> 

I The general concern of this paper is the dynamics of solid bodies in a fluid flow. This dynamics 

■ is relevant to many natural and industrial processes, like blood flows, sprays, or design of micro 

O i swimmers. 
■ 

^ ■ A main problem to understand this dynamics is to compute the drag exerted by the flow on the 

CN ■ bodies. From the mathematical point of view, a natural approach to this problem is to use the Euler or 

Navier-Stokes equations to model the flow. However, this generates serious difficulties. A famous one 
is D'Alembert's paradox, related to the Euler equation: in an incompressible and inviscid potential 
^ ■ flow, a solid body undergoes no drag |[2TI . 

I But the Navier-Stokes equations also raise modeling issues. Let us consider for instance a single 

solid moving in a viscous fluid. We denote by S{t) C M^, F{t) C M.^ the fluid and solid domains 
at time t, and := S{t) U F(t) the total domain. We assume that the fluid is governed by the 
Navier-Stokes equations. We denote up and pp its velocity and internal pressure, pp its density, i^p 
its viscosity. Thus: 

( pridtUF + UF -Vup) - HF^UF = -VpF - ppg, t > 0, X G 
I divnp = 0, t > 0, X G F{t), 
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with —ppg the gravitational force. In parallel to the fluid modeling, we write the conservation of 
lineal- and angular momentum for the body. Denoting xs{t) G the position of the center of mass, 
Us{t) G its velocity, and ujs{t) G the angular velocity at time t, these conservation laws read 



d f 
ms-j-Us{t) = - Y^Fvda - msg, 

"I JdS(t) 



— {Js^s(.t)) = -I {x -xs{t)) X {T.Fi^)da + PS {x -xs{t)) X {-g). 
dS{t) Js{t) 



dt 



Following standard notations, ps and ms := /9s|5(0)| are the density and mass of the solid (inde- 
pendent of t and x), x) E M3(R) is the newtonian tensor of the fluid: 

T,F = (2pfD{uf) - PF Id) , 
and Js{t) G M3(M) is the inertia matrix of the solid: 

Js{t) := PS / {\x - xs{t)\^ Id - {x- xs{t)) ®{x- xs{t))) dx. 
Js(t) 

The vector v = v{t, ■) is the unit normal vector pointing inside the solid S{t). Note that the velocity 
us{t, x) at each point x of the solid reads 

us{t,x) := Us{t) + ujsit) X (x-xsit)). 

To close the system, one usually imposes no-slip conditions, both at the fluid-solid interface and 
the cavity boundary: 

f upldsit) = us\ds{t) 
\ uplan = 0, 

and one specifies the initial data: the initial position of the solid Sq, 

UF,o ■= UF\t=o and ^5,0 := ^^5,0 + ^5(0) x (x-xso)- 

One could believe that system (ll.ll) - (ll.2b - (ll.3l ) is a good model for the interaction between a solid 
and a viscous fluid. Far from it: in the case of a sphere falling over a flat wall 

8(0) := 63+5(0,1/2), n := {^3 > 0}, 

it predicts that no collision is possible between the solid and the wall ! This no-collision paradox 
has been known from specialists since the 1960's, after articles by Cox and Brenner HI and Cooley 
and O'Neill Q in the context of Stokes equations. Since then, the no-collision paradox has been 
confirmed at the level of the Navier-Stokes equations (see |[T6l[r7l . and the preliminary result in Il23l ). 

Of course, such a result is unrealistic, as it goes against Archimedes' principle. It suggests that 
the Navier-Stokes equations are not relevant to collisional and post-collisional descriptions. Hence, 
many physicists have tried to find an explanation for the pai^adox. We shall focus here on one possible 
explanation, namely the no-slip condition. The idea is that, when the distance between the solids gets 
very small (below the micrometer), the no-slip condition is no longer accurate, and must be replaced 
by a Navier condition: 

f [uF - us) ■ v\as(t) = 0, [uF - us) X v\as(t) = -2Ps{^fv) x v\as{t), 

\ UF ■ v\dn =0, UF X v\dn = -2I3q {T.Ff) X z/|ao- 
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In other words, only the normal component of the relative velocity of the fluid is zero, to ensure 
impermeability. The tangential ones are non-zero, and proportional to the stress constraint, with 
constant slip lengths Ps-iPo. > 0. For a recent discussion of the Navier condition, notably in the 
context of microfluidics, we refer to |[T9l . See also the seminal paper |[T8l . Let us point out that the 
Navier-condition is sometimes used as a wall law, to describe the averaged effect of rough hydrophobic 
surfaces El . 

The effect of slip conditions (11.4b on collision was investigated recently by the authors in article 
|[T3l . More precisely, this article is devoted to a simplified model, in which 

• The Navier-Stokes equations are replaced by the steady Stokes ones (quasi-static regime). 

• The domain $7 is a half-space, the solid 5 is a sphere. 

In this context, denoting h{t) the distance between S{t) and the plane 90, it is shown that the dynam- 
ics obeys to the reduced ODE 

h = hV{h) + iPZ^g 
PS 

where the drag term satisfies = 0(| In as /i — )• 0. This is in sharp contrast with the 

case of no-slip conditions (11.31 ). for which ~ ^. In particular, it allows for collisions in finite 
time. We refer to |[T3l for all details and other results in the context of rough boundaries. 

Hence, paper |[T3l provides a resolution of the paradox: one can a priori keep the Navier-Stokes 
equations, up to consider the Navier boundary conditions (11.41) . Nevertheless, the analysis in |[T3l 
is limited to simple configurations and to Stokes flows. In the context of the full 3D Navier-Stokes 
system (11.11) . more complicated behaviors may occur. For instance, smooth solutions may exhibit 
singularities prior to any collision. To describe the qualitative features of the collision, one needs to 
consider weak solutions. The theory of weak solutions is well understood in the case of no-slip condi- 
tions and many references will be given in the next section. However, the existence of weak solutions 
with Navier conditions has been so far an open question, due to serious additional mathematical dif- 
ficulties. To address this question is the purpose of the present paper. Broadly, we shall build weak 
solutions for system (ll.ll) - (ll.2l) - (ll.4l ). up to collision between the solid and the cavity 0. 

The paper is organized as follows. Section |2] contains the statement of our main result: we give a 
definition of weak solutions, and state the existence of such solutions as long as no contact occurs. We 
explain the main difficulties in proving their existence, in comparison to the results available for no- 
slip conditions. We conclude Section [2] by an outline of our proof, to be carried out through sections 
[3] to [5] More precisely: 

• Section[3]is devoted to an auxiliary nonUnear transport equation, which is crucial to our approx- 
imation procedure. 

• Section |4] is dedicated to the construction of solutions for well-chosen approximations of the 
Navier-Stokes / solid dynamics. 

• Section |5] describes the limit procedure, from the approximate to the exact system. 



3 



2 Main result and ideas 



2.1 Weak solutions of Navier-Stokes with slip conditions 

The aim of this paragraph is to define a weak formulation and weak solutions for system (ll.ll )- (ll.2l )- 
(11.41 ). that is in the case of slip conditions of Navier type. We remind that in the case of no-slip 
conditions, the theory of weak solutions has been successfully achieved over the last ten years, first 
up to collision (see HI) and then globally in time (see ll23l in the 2D case, lITOl in the 3D case). Let us 
also mention the alternative approach in [3 ], and the recent result 1141 on the uniqueness of 2D weak 
solutions up to collision. 

As usual, in order to derive a weak formulation, the starting point is formal multiplication by 
appropriate test functions. These test functions must of course look like the solution itself. Notably, 
they must be rigid vector fields in the solid domain S. This forces the space of test functions to 
depend on the solution itself: it is a classical difficulty, already recognized in the no-slip case. A 
key feature of the slip conditions is that these test functions, and also the solution, will be moreover 
discontinuous across the fluid/solid interface. Indeed, the first Une of (11.41 ) ensures the continuity of 
the normal component, but the tangential ones may have a jump. It is a strong difference with regards 
to boundary conditions (11.3b . and it will generate many difficulties throughout the paper. 

We first introduce some notations for the classical spaces of solenoidal vector fields. Let O be a 
Lipschitz domain. We set 

V„{0) := W e T^iO), div ^ = 0}, V^{0) := {ip\o, ^ G V^iW")] , 
Ll{0) := theclosureof inL2(0), hI{0) := H\0)nLl{0), 

H^{0) := theclosureofP^(O) in iJi(C') 

We remind that elements u of L^(0) satisfy u ■ u = Oin H~^/^{dO). 

We also introduce the finite dimensional space of rigid vector fields in M^: 

TZ := {(fs, (ps{x) = V + u X X, for some V u e M^}. 

Finally, we define for any T > the space Tt of test functions over [0, T): 

Tt := {y^GC^([0,r];L2(5^)), there exists v?ir G P([0, T); P,(ll)), 995 G ^([0, T); 7^) 

such that ip{t,-) = ifpit, •) on F(t), ^p{t,■) = ^s{t, •) on S{t), for all t G [0,r]|. 

Let us point out once again that this space of test functions depends on the solution itself through the 
domains S{t) and F{t). Let us also notice that the constraint ip{t, •) G L'^i^) encodes in a weak form 
the continuity of the normal component at dS{t): 

(pF{t,-) ■ V = fs{t,-) ■ ^ at dS{t), \ft G [0,r). 
Multiplying (11.11 ) by ip £ Tt, integrating over F{t), and integrating by parts, we obtain (formally) 



I PFUF-(PF- [ PFUF-dt(pF- [ PF UF <^ Uf ■ V fF + [ 2fJ.FD{uF) ■ D{ipF) 
dt JF{t) JF{t) JF(t) JF{t) 

I {T^FV) ■ fF + i^Fl') ■VF+ PF{-g) ■ VF, 

Jan Jds{t) JF{t) 
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where the normal vector v at the right-hand side points outside 0, resp. inside S(t). Using (11.41 ): 

/ iJ^Fv)- ipF = \ {up ^ v) ■ {^F ^ v), 

Jan ^Pf^ Jan 

/ (J^Fl^) ■ = -T^r i{uF - us) X 1^) ■ {{lPf - ifs) X 1^) + / iJ^Fl^)-^S 

Jas{t) '^ps Jas{t) Jas(t) 



Eventually, one can use (11.21 ) to write differently the last integral: tedious but straightforward manip- 
ulations yield 

/ (Sfz^) • = -4: [ Psus ■ + [ Psus ■ dt^ps + [ Ps{-g) ■ ^s- 
Jas{t) M Js{t) Js{t) Js{t) 

Combining the previous identities and integrating from to T yields eventually 

f-T r f-T r r-T 



- / PFUF- dt(pF - / PSUS ■ dt(ps + / PFUF0UF ■■ V(pF 

Jo J F{t) Jo Js{t) Jo JF(i) 

+ 11 2fiFDiuF) : D{^f) + 1^1 
Jo Jfu) ^Pn Jo 



^{t) Jo JS{t) Jo JF{t) 

cT r rT r- 

/ {uf X u) ■ {ifF X v) 

^(t) "^Pn Jo Jan 

+ ^ f [ {{uF -us) -X v) ■ {{ipF - ifs) -X v) (2.1) 
'0 JdS(t) 



2/35 

/ / Pf{-9) ■^F+ / Psi-g) • V5 

Jo Jfu) Jo Js(t) 



+ / PFUF,0 ■ fF\t=0+ / PSUS,0 ■ <fs\t=0 

Jf{o) Js{o) 

Equation (12.11 ) is a global weak formulation of the momentum equations (11.11 ) and (11.21 ). taking the 
slip conditions (11.41 ) into account. Setting ip = uin the above formal computations yields that, for all 

te[0,T] : 

/ ]-pF\uF{t,-)\^ + [ lps\us{t,-)\^ + [ [ 2fiF\D{uF)\^ds 
JF(t) ^ Js(t) ^ Jo Jf{s) 



+ ^/ / I^^X^r + ^T^/ / \{UF-US)XU\' {1.1) 

^Pn Jo Jan ^Ps Jo Jds(t) 



< 



/ / Pf{-9)-uf + / / Ps{-9)-us+ / Pf\uf,o? + / Ps\us,o\ 
Jo JF{t) Jo Js{t) Jn\So JSo 



This goes together with the conservation of mass, that amounts to the transport of S by the rigid vector 
field n^. It reads 

dtxs + div (usxs) = inn, xsit,x) := ls{t)ix), 
or in a weak form: for all G V{[0, T),V{Ti)), 

^ [ dt^ - r [ us-V^ = [ ^\t=o. (2.3) 

Js{t) Jo Js{t) JSo 

Pondering on these formal manipulations, we can now introduce our definition of a weak solution on 
[0, T). We fix once for all the positive constants ps, pF, Pf, Ps-, Pn- 
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Definition 1 Let and Sq (Z VI two Lipschitz bounded domains o/ffi^. Let upfl G usfl G Tl- 

such that upfi ■ v = usfi • v on dSo . 

A weak solution of (ll.ll) - (ll.2l) - (ll.4l ) on [0,T) (associated to the initial data Sq, upfi and usfl) is a 
couple {S, u) satisfying 

• S{t) dO. is a bounded domain of for all t £ [0, T), such that 

Xs{t,x) ■.= lsit)ix)eL^i{0,T)xn). 

• u belongs to the space 

St := {nG L~(0, T; (0)), there exists up G L^{0,T; Hl{n)), us G L°°(0,r;7e) 
such that u{t, •) = upit, •) on F{t), u{t, •) = us(t, •) on S(t), for a. e. t £ [0, T]| , 



where F{t) ■.= n\ S{t)forall t G [0,r). 

• Equation (12.11 ) is satisfied for all if G It- 

• Equation (12.31 ) is satisfied for all ip G T^{[0, T);T>{^})). 

• Equation (12.21 ) is satisfied for almost every t G (0, T). 

Let us conclude this paragraph by a few comments on this definition of weak solutions: 

1. As xs G -^°°((0, T) X 0), the integrals over S{t) in (12.31 ) are integrable with respect to time: 
namely, 



1 1-> / dt'^ = I xsdt^ and th^ us ■ = / xsus ■ 
Js{t) Jn Js{t) Jn 

belong to L^(0, T). Actually, by the method of characteristics, as us G L°°(0,T;7^) (rigid 
velocity field), it is easily seen that 

Sit) = cPtASo) 

for an isometric propagator (pt^g which is Lipschitz continuous in time, smooth in space. It 
follows that all integrals in equation (12.11) make sense. For instance, as dS{t) is Lipschitz for 
all t and fields uf,us,(Pf,^s have at least L^ff^ regularity, the surface integral 

{{up - us) X n) ■ {{ifF - ifs) X n) 

dS(t) 

can be defined for almost every t in the trace sense. Moreover, it defines an element of (0, T). 
This can be seen through the change of variable x = (ptfiiy)'- the surface integral turns into 

J {t,4>t,o{y)) J^(^T{y)dy, 

dSo 

where 

]{t, x) := {{uF{t, x) - us{t, x)) X u) • {{ipF{t, x) - (ps{x)) X u) 

and where 

Jac,(y) = ||[V</>t,o(y)]-'^(y)||2det(V(/.t,o(y))(=l) 
is the tangential jacobian (see |[T5l Lemme 5.4.1] for details). This clearly defines an element 

of Li(o,r). 
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2. Equations (12.11 ) and (12.31 ) involve fields up, us, ^f, defined over Q and such that 

u = {1- xs)uF + Xsus, ip = (1- xs)v^F + XS'fS, 
However, a closer look at equations (12.11 ) and (12.31 ) shows that they only involve 

Xsus, xf(1,V)'Uf, as well as xs{'i-,dt)(ps and XF{l,dt,V)(fF- 
In particular, they only depend on u and cp, not on the choice of the extended fields up, us and 

^F,^S- 

3. The condition u G L°°(0, T; L^(ri)) imphes that 

up ■ v\Qs(t) = us ■ v\Qs(t) for a.e. t 
all terms being again understood in the trace sense. 

4. It is easy to see that equation (12. 31 ). that is the transport equation 

dtxs + div {xsus) = in P'([0, T) x Tl) 

can be written 

dtXs + div (X5?i) = in P'([0, T) x (2.4) 

and implies 

dtXF + diy {XFu) = Oin P'([0,r) x H), XF(t,x) = XF{t){x) (2.5) 



(remind that = $7 \ 5(t)). More generally, one can replace u by any v £ L°°(0, T; L^C^)) 
satisfying 

v{t, •) • i^bsw = • ^\ds{t) = us ■ ^\ds{t) for a.e. t 

where the last equality holds in the space H~^/^{dS{t)) (see ITll Theorem 3.2.2]). Note that 
equations (12.41 ) and (12.51) should be replaced by 

dtPs + div [psu] = 0, dtPf + div (pju) = 

in the case of inhomogeneous soHd and fluid, with variable density functions ps and pj. See 
lITOl in the case of no-slip conditions. Extension of the present work (on a single rigid and 
homogeneous solid in a homogeneous fluid) to more general configurations will be the matter 
of a forthcoming paper. 

5. Noticing that 

D{u{t, •)) = D{us{t, •)) = in S{t), D{^{t, •)) = D{^s{t, •)) = in S{t) 
it is tempting to write (12.11) under the condensed form 

- [ [ pu-dt(p+ [ [ pu(g)u: D{ip)+ [ [ 2pfD{u) : D{ip) 
Jo Jn Jo Jn Jo Jn 

= "boundary terms" (2.6) 
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where p := pfXf + PsXs< coupled with the global transport equation 

dtp + div [pu] = in O. (2.7) 

This kind of global formulation, reminiscent of the inhomogeneous Navier-Stokes equations, is 
used in the construction of weak solutions with Dirichlet boundary conditions: cf |[23l . However, 
it is not valid here: due to the discontinuity of the tangential components of u and ^p, neither 
difp nor D{u) and D{lp) belong to L^(0). For instance, 

dt^ = Xpdt^F + Xsdt^s + us ■v{ifF - Vs) ^ds 

where 5gs is the Dirac mass along the solid boundary dS. This is why we keep the formulation 
(12.11) . distinguishing between the solid and the fluid part. 

6. The definition of a weak solution that we consider can not be satisfactory after collision. Indeed, 
we do not specify any rebound law. Moreover, in the case of Dirichlet conditions at the fluid- 
solid interface, expUcit examples show that the analogue of our weak solution is not unique 

Eg. 

2.2 Main result 

Our result is the following 

Theorem 1 (Existence of weak solutions up to collision) 

Let $7 and Sq <^ Q. two C^'^ bounded domains ofM?. Let up^ G usfl G TZ such that 

uf,o • V = usfi • V on OSq. 

There exists T > and a weak solution of (ll.ll) - (ll.2l) - (ll.4l) on [0, T) (associated to the initial data 
So, upfi and us^). Moreover, such weak solution exists up to collision, that is 

S{t) (E n for all t e [0,T), and lim dist{S{t),dn) = 0. 

The rest of the paper will be devoted to the proof of the theorem. Briefly, there are two main 
difficulties compared to the case of Dirichlet conditions: 

• The lack of a unified formulation such as (12.61 ). 

• The lack of a uniform bound on solutions u. 

These difficulties appear both in the construction of approximate solutions, and in the convergence 
process. 

Indeed, the approximation of fluid-solid systems is usually adressed by relaxing the solid con- 
straint, through a penalization term. In this way, one is left with approximate systems that are close to 
density dependent Navier-Stokes equations. Roughly, they read 

dt{pnUn) + div{pnUn®Un) + ... = penalization ^^^^ 

dtPn + div {pnUn) = 0. 

In the case of slip conditions, in which a global formulation of type (I2.6I )- (I2.71 ) already holds, to build 
such approximation is quite natural. But in the case of Navier conditions, this is not easy. 
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Once an approximate sequence of solutions {pn,Un) has been derived, Diriciilet conditions allow 
for uniform bounds on n„. This simplifies a lot of convergence arguments, notably with regards 
to the transport equation 

dtPn + div {pnUn) = 

to which the classical DiPerna-Lions theory applies straightforwardly m. Also, it helps to provide 
strong convergence of Un in L^((0,T) x 0). In short, the lack of bound on dt{pnUn) (due to the 
penalization term) can be overcome by considering the fields Ps{t)un, where Psit) is the orthogonal 
projection in H^{Q) over the fields that are rigid in a 5-neighborhood of S{t). One can show that 
Ps{t)un has good equicontinuity properties uniformly in S and n. 

In the case of Navier boundary conditions, no uniform bound is available in . This forces us to 
use more the structure of the solution n, in particular- the bounds on the fluid and solid domains 
separately. This is also a source of trouble for the construction of approximate solutions, as one must 
find an approximation scheme in which such structure is not too much broken. 

2.3 Strategy of proof 

Let us describe here briefly the main lines of our proof. Let So,uf,o,us,o as in Theorem [T] and 

PO ■= /Of(1-15o) + Ps'i^So, Uo := {l-lso)uF,0 + l5o'"S,0- 
The keypoint is to consider approximate problems of the following type: find {S"',u"') such that 

a) S'^{t) dVtis a bounded Lipschitz domain for all t G [0, T], such that 

xUt,x) := lsn(i)(x)GL-((0,T)xJ7)nC([0,r];Lf(J7)), Vp<+oo 

b) G L-(0, T; (J^)) n ^^(O, T; Hl{n)). 

c) For all if e H^{0,T; Ll{n)) n L'^{0,T; Hl{n)) s.t ip\t=T = O: 

^ /" p" {u^dtip + t;" ® u'' : V(^) + r [ 2^"Z)(n") : Diip) 
Jn Jo Jn 

I rT r ^ rT 



+ / Kxz.).(^xz.) + — / / (K-P>")xz.).((9.-P5»x^) 

^Pn Jo Jan ^Ps Jo Jas^it) 

+ n r [ xSK - • - PsV) = r [ • ^ + [ Pouo- ^\t=o 

Jo Jn Jo Jn Jn 

d) dtx^s + • VxS = 0, Xs\t=o = Iso- 

In above lines, 

• /?" := pf{1 — Xs) + is the total density function. 

• := pf{1 — Xs) + ^Xs is an inhomogeneous viscosity coefficient. 

• Pg = Pgit) is the orthogonal projection in L^(S'"(t)) over rigid fields. This means that: 
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VO<t<r, V^i5e7^, VnGL^(fi), P^{t)u£n and 



/ 



xWr){n-PS{t)u)-us = 0. 



• Eventually, 



is a field that satisfies 



u"(t, •) outside a 6 neighborhood of 5'"(t), t G [0, T) 



for some 5 fixed and arbitrary in (0,dist(5o, dQ)/2). Moreover, will be chosen so that it is 
close to u"^ outside 5" (in topology). In this way, it will asymptotically coincide with the 
limit M of u". Further details on the definition of i;" will be provided in due course. 

Let us make a few comments on such approximate problems: 

1. They rely on the use of the fields Pg U^, that were already introduced in [Si in the context of 
Dirichlet conditions. These fields appear both: 

i) in the transport equation for Xs- They will allow for a good control of the trajectories of the 
approximate solid bodies 5"". 

ii) in the penalization term n / Xsi'^"^ ~ -Ps^") • ~ ^s^)- Formally, as n goes to 



infinity, this term will allow to recover the rigid constraint inside the solid. 

2. Note that a contrario to the large penalization term, the viscosity term vanishes asymptot- 
ically in the solid part. Hence, there will be no uniform bound in H^{^) for u", as expected 
(see the discussion in paragraph 12.21) . 

3. A specificity of these approximate problems is that the transport equation d) is nonlinear in 

for a given u^. Indeed, Pg depends on X5 (c/the formula in section O. The whole sec- 
tion [3] is dedicated to this auxiliary nonlinear transport equation, which is a keystone of the 
approximation procedure. 

4. Once the solution of d) is found and seen as a functional of u^, equation c) can be written 
as F{u^) = for some functional F from L°°L^ n L^H^ into itself. In short, we shall solve 
this equation by a Galerkin procedure: we shall look for an approximate solution u"'^(t, x) = 
Ylk=o ^k{t)ek{x) where (cfc) is an orthonormal basis of L^(r2). We shall solve approximate 
equations J^'^ (u^'^) = by Schauder's theorem and pass to the limit with respect to N. This 
process is explained in section ID 

5. Note that the field f " satisfies 



In particular, one can write 



dtx's + ■ VxS = 0, and Otp" + • Vp" = 0. 
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This will allow to obtain energy estimates in a standard way, in the spirit of the approximate 
systems (12.81) used for Dirichlet conditions. The price to pay is the necessary control of — v", 
which will exhibit strong gradients near dS" as n — > +00. Moreover, a similar "boundary layer 
behaviour" will be involved in the approximation of discontinuous test functions (/? G 7t by 
continuous test functions if"' (involved in c)). The whole convergence process will be analyzed 
in section [21 

3 A nonlinear transport equation 

Let T > 0, n G L°°(0, T; LK^I)). This section is devoted to the equation 

dtXs + Psu • Vxs, Xs\t=o = Iso, 
where Psu is defined by the following formula 

Psu := j^PsXsu + (^J^^ J^psxs {{x - xs) X u) dx')^ x {x - xs) (3.1) 

where the center of mass, total mass and inertia tensor of the solid are defined by 

xs ■= / PSXS, M := PSXS, (3.2) 



If xs{'t,x) = ls(t){x) with S{t) a subdomain of Q, , Ps{t) is the orthogonal projection in L'^{S{t)) 
over rigid vector fields, see IS. 

We start with the regular case, that is when u G C([0, T]; P^l^))- This case will be useful for 
Galerkin approximations of a)-d). 

Proposition 2 (Well-posedness) 

Letue C{[0,T];V^{Tl)). 

i) There is a unique solution xs e L°°{{0, T) x R^) n C([0, T];LP{^^)) {p < 00) of 

dtXs + div {xs Psu) = in xs\t=o = Iso- (3-4) 

//) Moreover xsit, ') = ^ s{t) for all t, with S{t) a Lipschitz bounded domain. More precisely, 

Sit) = <Pt,oiSo) 

for the isometric propagator cj)t^s associated to Psu : [t, s) 1— )■ <f)t,s G C'^ClO, T]'^; Q'^(R^)). 
Proof. We can suppose u G C([0, T];V„{R^)) with no loss of generaUty. 



Assume for a moment that we have found a solution xs of (13.41 ). Then, we can see (13.41 ) as a linear 
transport equation, with given transport Psu G C([0,T]; TZ). By the method of characteristics, we 
get easily 

xs(.t,(f>t,o{y)) = '^soiy), (3-5) 
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where (f)t,.s is the isometric propagator defined by 



(3.6) 

dt(i)tAy) = Psu{t,<Pt,s{y)) v(s,t,x) G (o,r)2 xm3. 



Now, we use (13.51) in the expression (13.11) for Psu. We obtain: 



Psu{t,x) = ^ [ psln{(t>tfl{y))u{t,(t>tfi{y))dy 
-'So 



+ (^J-^(t) ^ P5 in(<At,o(y)) (<At,o(y) - xs{t)) X (/.t,o(y))t^y) x - xs{t)) (3.7) 

where M := \So\ps, xs{t) := / ps(j)t,oiy)dy, and 

^(i) := ^ P5(|<^t,o(y)-a;5(t)P/d-('/'t,o(y)-2;5(t))®('/'t,o(y)-X5(t)))dy. 

In particular, denoting Isom(M^) M"^ x 03(M) the finite dimensional manifold of affine isometrics, 
we deduce from (13.61 ) and (13.71 ) that t i— )• (ptfi, [0,T] i— )• Isom(M^) satisfies an ordinary differential 
equation, of the type 

^'/>t,o = Us{t, 4)tfi), 00,0 = Id, (3.8) 

for a time-dependent vector field Us over Isom(M^). Namely, Us{t,(j)) G T(^(Isom(M'^)) w 7^ is 
defined by the same formula as in (I3.7I ). replacing everywhere o by 

Conversely, if we manage to show existence of and uniqueness of a solution of (13.81) over 
[0, T], then formula (13.5b will define the unique solution xs of the nonlinear- equation (13.41 ). proving 
Theorem 121 

Hence, it only remains to study the well-posedness of (13.81 ). We can identify Isom(M'^) with 
X (93 (M) C X M^, and identify all tangent spaces with 7^ C x M^. By the Cauchy-Lipschitz 
theorem, there is existence and uniqueness of a maximal solution if Us is continuous in t, cj), 
locally Lipschitz in 0. Considering the expression of Us, see (13. 71 ). this follows from 



Lemma 3 Letv G C{[0,T];C^^{R^)). Then, the function 

M : [0, T] X Isom(]R3) ^ R, 0) = / lf^(</,(y))i;(t, 



5" 

w continuous in {t, (j)), and uniformly Lipschitz in <j) over [0, T]. 

Proof of the lemma. The continuity is obvious. Then, for two affine isometrics cj) and (f)' , we write 

M{tA)-M{t,ck')= f Inim) {v{t,cl))-v{tA'{y))) + I {iniHy)) - '^n{c^\y)))<tA'{y)) dy 
J So J So 

:= Mi{t)+M2{t). 

Clearly, 

|Mi(t)| < sup \dMt,x)\ [ \cP{y)-(l)'{y)\dy < C^,^'\\<p - (p'Woo- 
te[o,T], .750 

l^l<li(0,0')l|oo 
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As regards Af2, we write 

M2(t) < sup \v{t,x)\ I \ln{c^iy))-lnicl)'{y))\dy < C^, [ \lMy)) - lni(l>' {y))\ dy 
ie[o.r], 

kl<ll0'l|oo 

For each y, the integrand is non-zero if and only if (f){y) G Q and (/>'(?/) G Q'^ or vice-versa. As 
\4>{y) — cl)'{y)\ < 1 — 0'||oo, this is only possible if cl){y) and (p'iu) are in a ||0 — i;^'||oo -neighborhood 
(say V) of dQ. Hence, 

\M2{t)\ < 

This concludes the proof of the lemma. 

Last step is to prove that the maximal solution is defined over the whole interval [0, T]. From 
(I3.6I )- (I3.7I ). one can write 

(ptfiiy) = xs{t) + Qs{t)y 

where xs{t) is defined in (13.61) and Qs{t) is an orthogonal matrix. In particular, the only way that the 
maximal solution is not global on [0, T] is through a blow-up of xs- But, again, from (13.71) . 





dy] 







. d , . . 1 



PS ^Q{(t>t,o{y)) u{t, (f)t^o{y))dy 

So 



< C ||n||£oo((o,T)xn) 



which prevents any blow-up. This ends the proof of the theorem. 
Proposition 4 (Strong sequential continuity) Assume that 

vJ" in C([0,r];P^(17)). 

Then with obvious notations, one has 

X's ^ XS weakly * in L°°((0, T) x M^), strongly in C([0, T]-Ll^{W')) (p < oo), 
as well as 

P^n" ^ Psu strongly in C([0, T]; Q~ (M^)), 0- ^ <^ strongly in C\[^, r]^; Q~ (M^)). 

Proof of the proposition. As n" converges in C([0, T]; 'Dcr(^)), we have that PgvT' is bounded in 
L°°(0, T; L^(r2)) . Similarly is bounded in L°°((0,r) x Q). Furthermore, up to a subsequence 
that we do not relabel, Pgu"' converges weakly-* in L~(0, T; ^/^^(M^)) to some us and Xs(0)(= 
l^p, for all n G N) converges strongly in L^{Q.). Applying Di Perna-Lions theory, we obtain that Xs 
converges weakly-* in L°°((0, T) x Q) and strongly in C([0, T]; Lf^^($7)) for all finite p. Its Umit xs< 
satisfies : 

dtxs + div{xsus) = 0. 



Using the convergence of both Xs ^^'^ i'^ equation (13. IK we obtain that us = Psu, where 
Ps is defined similarly to Psu, replacing xs by Xs- Moreover, the convergence of P'gU holds in 
C([0, T]; C;'^(M^)). Consequently, {xs,Psu) is the unique solution of (13.41 ) so that xs = Xs and 
Psu = Psu, and all the sequence converges. 

To derive the convergence of the propagators from the convergence of the vector fields PgU^ 
is then standard, and we omit it for brevity. 
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Proposition 5 (Weak sequential continuity) 

Let (n", X5) be a bounded sequence in L°°(0, T; L^(rj)) x L°°((0, T) x VL), satisfying 

r/ie«, up to a subsequence, one has 

u"" weakly * in L~(0, T; LI{^)) 

Xs^XS weakly * /« ( (0 , T) x ) , strongly in C{[0,T]; Lf^^{R^)) (p < 00), 
with (us, Xs) <^ solution of 

dtXs + div (Psuxs) = inR^, xs\t=o = Iso- 

Moreover, xs satisfies condition ii) of Proposition^ and the following additional convergences hold: 

P^n" ^ Psu weakly * in L°° (0, T; C^^ (M^) ) , 

(/)" ^ weakly * in T^^'°°((0, Tf-C^^{R^)) strongly in C([0, Tf- C^^{R^)) ■ 

Proof. The proof follows the same scheme as the previous one. We only sketch the arguments. First, 
up to the extraction of a subsequence, we obtain that 

u"" weakly * in (0, T; (f]) ) 

Then, as before, we obtain that P^n" is bounded in L~(0, T; (E^)). This yields that 

P^u" ^ Pu weakly * in L°°(0, T; F^ocll^^)) 

(still up to a subsequence). We then deduce applying Di Perna-Lions theory that, up to the extraction 
of a subsequence, X5 converges strongly in C{\Q,T]] L^^^iR^)) to some xs> which in turn implies 
that Us = Psu and that {xs, Psu) is a solution to our tranport equation. Eventually, uniform bounds 
on and (which imply weak-* convergence of a subsequence in W^''^) follow easily. 

4 Approximation 

This section is devoted to the resolution of approximate fluid-solid systems. These approximate sys- 
tems were introduced in paragraph 12.31 c/a)-d). The previous section has focused on the transport 
equation d). It remains to examine c). At first, we explain a little how the field connecting P^u"' 
to is defined. The detailed definition of f " will be achieved in section |5] 

4.1 Connecting velocity 

We first remind a classical result on the equation div u = f, taken from ifTTl Exercise III.3.5]: 

Proposition 6 Let O be a bounded Lipschitz domain. Let f G and ip G H^/'^{dO) satisfying 

the compatibility condition 
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Then there exists a solution u G H^{0) of 

div u = f in O, u = ip at dO 

with 

\\'^\\m{o) < Co {\W\\m/-i{do) + 11/1^2(0)) ■ 
The previous proposition yields easily 

Corollary 7 (Extension of solenoidal vector fields) 

There exists a continuous linear operator Eq : H^{^) i— )• H^{M.^) satisfying Equ = u on Q. 
Moreover, for all open subset a; ^ fi, 

ll^!^'"ll/fi(R3\aj) < C't^ll'"ll/fi(n\u)! y u e H^i^l). 

Corollary 8 (Connection of solenoidal vector fields) 

For all 6 > 0, there exists a continuous linear operator 

V' : H^, \ So) X H',(S~o) ^ H^R^), {U, Us) ^ V'[U, Us] 

such that 

V^[U,Us] = Us in So, 

V^[U,Us] = U outside a 5 neighborhood of So- 
From there, we have the following 

Proposition 9 For all 6 > 0, there exists a continuous mapping 

: L\0,T;Hl{W')) X L°^{0,T;TZ) ^ L\0,T; hI{R^)), {u,us) ^ v^[u,us] 

such that 

v^[u,us]{t,-) = us{t,-) in 5(i), 

v^[u, us](t, •) = u{t, •) outside a 5 neighborhood of S{t), t € [0, T), 

where, as usual, S{t) := (ptfiiSo) and (j) = (j)t,s the isometric propagator associated to us- 
Moreover, can be chosen so that 

\\v^[u,Us]\\l2(^0,T;HHR^)) < ^ I ') II J^i (k3\5(^) + \Ws{t, ^Wli^sit))) dt, 



where C depends on 5 and T. 

Proof of the proposition. The proposition can be deduced from Corollary [8] using Lagrangian 
coordinates. Namely, we introduce U and U s through the relations 

u{t,(ptfi{y)) = d(j)tfi\y{U{t,y)), us {t,(l)t,o{y)) = d(ptfi\y {Us{t,y)) . 
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Clearly, for all t, U{t, •) and Us{t, ■) define elements of i?^(M^ \ ^o) and H^{So) respectively. Using 
Corollary [8j we define v^[u, us] through the relation 

v'[u,us]{t,My)) = d(t>t,o\y{v'[Uitr),Us{t,-)]{y)) ■ 
It fulfills all requirements, which ends the proof. 
Back to system c), the idea is to define 

:= v^[Enu'',PSu'']. 

Clearly, for any time T" such that 

dist(5"(t),ai7) > 25, tG[o,r"], 

u"|n will belong to L'^{0, T"; H^{n)) and will satisfy 

v^it,.) = PSit)u^it,-) in 

.) = •) outside a 5 neighborhood of t G [0, T"). 

Let us stress that there is still some latitude left in the construction of t;", through the choice of the 
operator in Corollary [H As will be shown in section [51 this operator can be chosen depending on n 
(V^ = y^.") so that i;" is close to n" outside 5" (in topology). However, this additional property 
will not be needed until section Hi 

Last remark: the resolution of a)-d), and the whole construction of weak solutions, will be first 
performed on a small time interval [0, T], for a time T that is uniform in n. Existence of weak solutions 
up to collision will follow from a continuation ai^gument, to be explained at the end of section |5] 

4.2 Galerkin approximation 

As pointed out in paragraph 12.31 the resolution of a)-d) is canied out through a Galerkin scheme. Let 
(efc)fc>i being both an orthonormal basis of L'^{0,) and a basis of H^{Q), with elements in Va-i^)- 
The aim of this paragraph is to find for all N, n and some T > a couple (5"^, n^) satisfying 

a') (t) C O is a bounded Lipschitz domain for all t G [0, T], such that 

X^{t,x) := l5iV(,)(x)GL-((0,r)xf^)nC([0,r];LP(17))(p<oo) 

N 

b') u^{t,-) = ^afc(t)efc, with a = {ai,...,aN)eC{[0,T]f. 

i=l 

c') For all ip £ V{[0,T); span(ei, . . . ,eAr)) 

^ / (n^ • dtip + v''®u^ : Vv9) + T [ 2ii''D{u'') : D{ip) 
Jn Jo Jn 



+ ^/ / xz.).((^xz.) + — / / ((u--P^n-)x^).((v.-P^9')xz.) 

^Pn Jo Jan ^Ps Jo Jds^(t) 

+ n r [ (tx^ - P#n^) • (ip - P^ip) = r [ p^'i-g) ■ + [ pouo- ip\t=o 

Jo Jn Jo Jn Jn 
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d') dtx^ + p^tx^ • vxf = in n, \t=o = iso- 

In above lines, similarly to the original problem: 

• := pf{1 — Xs) + PsXs is the total density function. 

• /i^ := — Xs) + '^Xs is an inhomogeneous viscosity coefficient. 

• Pg = Pg (t) is defined by (13. IK adding the upperscript everywhere. 

• Eventually, v'^ = v^[u^ , P^ u^], see paragraph 14. 1 1 

Note that all quantities above depend on n, notably through the penalization term and the viscosity 
coefficient. But we omit n from the notations to lighten writings. Also, note that can be seen as 
an element of L?{Q, T; H^{M^^)), as the are defined globally. In particular, = v^lu'^ ,Pg u^] is 
well-defined. 

The main result of this paragraph is 

Theorem 10 There is T > 0, R > 0, such that for all n,N, a')-d') has at least one solution such 
that \\u^\\L°°(o,T;L^n)) < R- 

To prove Theorem [TOl we shall express our Galerkin problem as a fixed point problem, and will apply 
Schauder's theorem to it. Thus, we want to identify as the fixed point of an application 

■T-N 

T : ui-?- u, 

defined on := {u & C{[0,T]; span{ei, . . . ,eN)), ||ii||^cx)(o^T;L2(Q)) < -R}. We proceed as 

follows. Let u G Br^t- 

• Step 1. Let be the solution of 

dtXs + Psu ■ Vxs = 0, xs\t=o = Iso, 

given by Proposition |2] We know that xsi^ix) = ls[t){x) with S{t) a bounded Lipschitz 
domain, t G [0, T]. We define accordingly: 

p := pf{1-Xs) + PSXS, P ■■= fJ'F{l-Xs) + -^XS, v{t,x) := v^[u,Psu]. 

• Step 2. We consider the following ODE, with unknown u : [0, T] i— )• span(ei, . . . ,6^)- 

^(t)^n(t) + B{t)u{t) = fit), n(0) = := ^ ( / no • efc ) e^, (4.1) 

k=i ^-^^ ^ 

in which A{t) := {aij(t))i<ij<N, B{t) := {bi,j{t))i<ij<N and f{t) := {fi{t))i<i<N are 
defined by 

pCi • 6j , 

hj '■= I P{v ■ ^^j) -e-i + I 2/xZ:>(ei) : D{ej) + / (e^ x u) ■ {ej x u) 
Jn Jn ^Pn Jan 

+ ^ / ((^i ~ ^sGi) X v) ■ {{ej - PsCj) X u) + n I xsi^i - Ps^i) ■ (ej - Psef 
'^Ps Jds{t) Jn 

fi ■= / pi-a) ■ ei. 
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We have identified liere the function u with its coefficients in the basis ei, . . . , bat. Note that 
the function p defined in step 1 has a positive lower bound, so that A{t) > mm{ps, Pf)In as 
a symmetric matrix, whatever the value of xs- Also, the continuity of A and B over [0, T] is 
easy and will be proved below. In particular, equation (14.11 ) has a unique solution 

u£ CH[0,r]; span(ei,...,e7v)). 

In this way, we can associate to each u € Br^t some field 

u = F^{u) e C([0,T]; span(ei,...,eAr)). 

The whole point is to prove 

Proposition 11 There exists T > 0, R > 0, uniform in n and N, such that is a well-defined 
mapping from Bh^t to itself continuous and compact. 

Before proving this proposition, let us show how it implies Theorem [TO] By Schauder's theorem, it 
yields the existence of a fixed point G Br^t of . Let = I^jv be the corresponding solution 
of the transport equation on [0, T] x M'^. As will be clear from the proof, the time T of the proposition 
satisfies 

dist(5^(t),aO) > 25, VtG[0,T], 

for some 5 fixed and arbitrary in (0, dist(S'o, dVL) /2). Hence, a') is satisfied, and := [u^, Pg u^] 
satisfies ■ u\q^ = {), d& well as 

dtp"" + • Vp^ = in 

(see remark 4 after the definition of weak solutions, and remark 5, pai'agraph 12.31 ). Finally, we notice 
that ODE (l4~n ) is equivalent to: for all ip G P([0, T); span(ei, . . . , cn)) 



''^ [ p^'dtu^'-v + 



/o Jn 



[ I p^t;^ . Vn^ • + /" I 2^^Z)(n^) : D{^) 
JO Jn Jo Jn 



+ / K xz.).(v.xz.) + — / / ((n^^-Pi^n^^)x^).((^-P^V,)x^) 

2pn Jo Jdn ^Ps Jo Jds'^it) 

Jo Jn Jo Jn Jn 

(4.2) 

Combining this equation with the previous one on p^ leads to c')- Note that condition ■ uIqq = 
is needed for the convective term to vanish through integration by parts. 

Proof of the proposition. 

Step 1: Definition of J-^ . 

We first prove that is well-defined from Bji t to C([0, T]; span(ei, . . . , cat)) for any T and 
it! > 0. The only thing to check is the continuity of matrices A and B in (14.11 ) with respect to 
time, which will guarantee the existence of a solution to the linear ODE (14.11 ). As xs belongs to 
C{[0,T]; LP{Q)) for all finite p, so does p, and A is clearly continuous. As regai^ds B, the only 
difficult terms are 

lit) := / pivit, •) • Ve,) • e„ J(t) := 77^ / ((e^ - Psit)ei) x z.) • ((e,- - Psit)ej) x i/). 
Jn ^Ps JdS{t) 
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We remind that the propagator (f) = (j)t^s associated to Psu satisfies 

</.GCi([o,r]2;Q-(]R3)) 

Hence, a look at the construction of , c/ Corollary [8] and Proposition |9] (see also Lemma [T6] in the 
appendix IaI). yields 

7;GC([0,r];Fi(]R3)). 

It implies that 1 1— I{t) is continuous. 

As regards J{t), we change variables to go back to a fixed domain. We set x = (ptfiiu) to obtain 

^Ps Jaso 

where 

:= {{ei{x) - Ps{t)ei{x)) x u) ■ {{ej{x) - Ps{t)ej{x)) x u) 

and where 

Jac^(y) = \\d<pt,o\y^ J^{y)\\2det{d(l)t,o\y){= 1) 

is the tangential jacobian. See 03 Lemme 5.4.1] for details. As j is continuous in t and smooth in x, 
we obtain that 1 1— )■ J{t) is continuous. 

Step 2: sends Br^t to itself. 

Here, we need to restrict to small T. More precisely, we fix < < ^dist(S'o, d^), and consider 
a time T such that 

inf dist(5(t),(90) > > (4.3) 

Let us prove that such time T does exist and can be chosen uniformly with respect to N and n. For 
all u G Br^t, we write 

S{t) = (ptfl{So) 

with (f) the propagator associated to the rigid field Psu = X5 + a;5 x — xs) defined in (13.11 ). It is 
enough that 

,^ , , dist(5o,90) - , 
sup \dtct>tfl{t,y)\ < — ^ ^ , t e [0,r], y e 5o. 

tG[0,T] J- 

Wefind 

\dt(t>tfi{t,y)\ < \us{t,(l)tAt,y))\ <\xs{t)\ + \ujs{t)\\y - xso\ 
using that the propagator is isometric. Moreover, classical calculations yield 



\xs{t)\^ + J{t)usit)-u;sit)= [ ps\Psuit,-)\^ < [ Ps\uit,-)\^ < PsR^- 

Js{t) Js{t) 

We can then use the inequality 



\xs{t)\ + \ujsit)\\y - xso\ < V2 max{l,\y - xs,\) {\xs{t)\^ + \uJs{t)\^) 



< 



Co {\xs{t)\^ + J{t)ujsit)-ivs{t)y^^ 
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where for instance 

max(l,sup gs -x^J) 

c7o := v2 ^ — , ^, Ao : smallest eigenvalue of J(0). 

min(l, Ao) / 

dist(5o,517) -2(5 .„ . , 
Eventually, any T < ^ ^ ^'.^ will satisfy (|43] ). 

Let now ti be arbitrary in Thanks to (14.31) . we have that v = v^[u,Psu] satisfies u-z^I^q = 0, 

and 

dtp + V ■ Vp = in il. 

Mutiplying (14.11 ) by u, integrating in time, and combining with the last transport equation, we obtain 
the energy estimate 

\\^pHt,-)\\l, + f f 2p\D{u)\'' 
Jo Jn 

+ TiT / + / \{u - Psu) X + n Xs\u-Psu\^ (4.4) 

^Pn Jo Jan ^Ps Jo JdS(t) Jo Jn 

< [ [ p{-g)-u + [ po\u^\^ 
Jo Jn Jn 

As min(pi;', ps) < p < m.ax{pF, ps), we deduce easily that 

II^^IIl°°{0,T;L2(Q)) < R 

for R = R{T, uq) lai^ge enough. Hence, sends Bji^t to itself. 

Step 3. Compactness of . 
N 

For any u = a^e^, we get from equation (14.11 ): 

k=l 

\ja{t)\ < \A-\t)\\B{t)\\a{t)\ + \f{t)\ < R\A-\t)\\B{t)\ + \f{t)\. 
Integrating with respect to time, we obtain 

sup (\a{t)\ + \^a{t)\] < C 
te[o,T] V at J 

(where the constant at the r.h.s. may depend on N or n). In other words, 

sup ||-7^^(n)||cl([0,T];Span(ei,...,ejv)) < C" 

which provides compactness in Bji^t by Ascoli's theorem. 
Step 4. Continuity of . 

Let {u^) a sequence in Br^t, such that — )• u in Br^t (that is uniformly over [0,T]). We 
want to show that {u^) — F^{u) in Br^t- First, we note that, as span{ei, . . . , cat) is a finite- 
dimensional subspace of Va-{0.) we have that u'^ converges to u in C{[0,T];V^{ri)). Then, we use 
Proposition m With obvious notations, 

xl ^ xs weakly * in L°°((0, T) x W"), strongly in C{[0, T];LPJW')) {p < oo), 
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as well as 

pIu'' Psu strongly in L°°(0, T; C^^{W^)), ^ 4> strongly in VF1'°°((0, Tf; (M^)). 

From there, and the construction of (Corollary [8j Proposition |9j Lemma [TT] in appendix lAl). it is 
easy to see that 

strongly in C([0, T]; //^(R^)). 



By slightly adapting the arguments of Step 1, one can then show that the matrices in (14.11) satisfy 

^ B, ^ A strongly in C([0, T]). 
From classical results for ODE's, it follows that 

= J^{u'') u = F{u) strongly in C([0,r]; span(ei, . . . ,eAr)). 
For the sake of brevity, we leave the details to the reader. 

4.3 Convergence of the Galerkin scheme 

In the previous paragraph, we have built for each n,N a solution u^'^ (denoted for brevity) of 
a')-d')- It is defined on [0, T] for some time T uniform in n, N, satisfying (14.31) . The next step is to let 
N go to infinity, to recover a solution of a)-d). We remind the uniform energy estimate (see (14.21) ) 



^ ' Jo Jn ^Pn Jo 
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+ ^ [ [ \{u^- P^u^) X + n / [x§\u^- Ps^u^l' (4.5) 
^ps Jo Jds'^it) Jo Jn 

< f f p{-9)-u^ + [ Po\u^\' 
Jo Jn Jn 

It yields that 

{u^) NeN is bounded uniformly with respect to N in L°°(0, T; Ll{n)) D L^{0, T; Hl{n)) 

The bound in follows from the bound on L'(n^) and Korn's inequality, see Il22l . From there, 
we will be able to show strong convergence both in the transport equation d') and in the momentum 
equation c')- As regards the transport equation, we rely on Proposition |5] Up to a subsequence, one 
has 

weakly */weakly in L°°(0, T; L^(J])) n L^{0, T; Hl{n)) 
for some u{= u"), and it follows from this proposition that 

Xs ^ XS weakly * in L-((0, T) x R^), strongly in C{[0, T];Ll^(R^)) (p < oo), 

as well as 

P^u^ Psu weakly * in L~(0, T; Cj^^iM.^)), 

(p^^cp weakly* in VF^'~((0, r)^; C7~ (M^)) , strongly in ^([0, T]; C;'^^(M3)) . 
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up to another extraction. We stress again that all limits depend on n. 

It remains to study the convergence of equation c')- Therefore, we fix the test function: we take 

ip{t,x):=x{t)ej, xG^([0,T)). 

for some fixed j. The point is to obtain as — >■ +00 the limit equation c), still with ip{t, x) = 
x{t)ej{x). But as j is arbitrary, and as {ek)k>i is a basis of H^{Q,), standard density arguments will 
allow to extend the formulation to general test functions. 

At first, we need to prove that, 

v^'ln ^ v = v^[Enu,Psu]\n in L\0,T- H^in)). 
It is enough to prove that 

t)^ := /[Enn^,P^tx^] ^ V := v^[Enu,Psu] 

weakly in L^(0, T; Hi^^{M.^)). In view of Corollary [8] and Proposition |9j it is an easy consequence of 
Lemma [TTl 

We are now ready to handle the asymptotics of c') (with ip{t, x) = xit) ej{x)). As before, for the 
sake of brevity, we focus on the two most difficult terms, those which involve 

Jn 



As regards J^(t), once again we change variables to go back to a fixed domain. We obtain 

1 

2^ 



dSo 
where 

f{t,x) := - X u) ■ ((e,(x) - P#(t)e,(x)) x 

and where 

Jacf(y) = ||d<ol,"'Ky)ll2det(d<Ai^ol.) = l- 

Let := Equ^ — Pg , resp. ryj^ := ej — Pg ej to which we associate R^, resp. Hj^ through 
the change of coordinates: 

r^{t,cp^{t,y)) ■.= d<p^\yR^{t,y), rj^ (t, it,y)) := d<pf\yH^ {t,y) . 

From the weak convergence of , we deduce that converges weakly in L'^{0,T; Hj^^{M.^)). 
Given the strong convergence of in C([0, T]; we also have that ryj^ converges strongly to 

rjj := ej — Ps&j in L^(0, T; //^^^^(R^)). Furthermore, as d(f)fQ\y is an isometric mapping for all N, 
we get that : 

f{t,ct>tfl[y)) = (i?^ X y) ■ {Hf xu), ViV G N 
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where, because of lemma [TT]: 

R'^ ^ R weakly in ^^(0, T; H'^{n)) , Hf Hj strongly in ^^(O, T; H^{Q)) . 

with obvious notations. This yields corresponding weak and strong convergences of the traces of these 
functions on OSq . Having in mind that Jac^ = 1 for all N, and going back to the moving domain, we 
obtain easily that J converges weakly inL^(0,r) to : 

1 



J{t) :-- 



dS{t) 



{{u - Ps{t)u) X u) ■ {{ej - Ps{t)ej) x u). 



We finally turn to the convergence of , for which we will need some compactness on (p^u^). 
Therefore, we introduce some notations: we denote by P the orthogonal projection from onto 
L^(r2), respectively Pk the orthogonal projection from L'^{Q) onto span(ei, . . . , Cfc). We also remind 
that our strong, resp. weak, convergence results on p^, resp. imply that. 

p^u^ pu weakly-* in L°°{Q,T]L^{Q)). 

In particular, we have for any fixed k: 

Pk{p^u^) Pk{pu) weakly-* in L°°(0,T; Ll{n)) as N 

Moreover, equation c') can be written: for all 1 < A: < A^, 

atPfc(pV) + PfcF^ = in V'{Q,T-[Hlm*) 

where G V (O, T; [HI{Vl)]*)) is defined by the duality relation: 



oo. 



(4.6) 



Jn 



+ 



1 



T 



2/35 Jo Jds^{t) 



+ n 



Jo. 



x^{u^-P^'nn-{y^-P^\) + 



:Vip - [ [ 2p.^D{u^) : D{^) 
Jo Jn 

((n^-P#^.^)x^).((c^-P|^^)x^) 

1 



N, 



{u^ X I/) ■ {ip X 1^) 



Jn 



2/3n Jo Jan 
foTa\l^eV{0,T;H^{n)) . 



We remind that for / G [hI{Vl)]* , Pk is defined by duality: < Pkf.'P > := < f.PkV >■ From 
the above expression for and the various bounds already obtained, it is easily seen that for any 
fixed k, (PfcF^) is bounded (in N) in L^{0,T; [H^in)]*). Hence, the same conclusion applies to 
{dt Pk{p^u^))- Combining with (I4.6I ). it follows that for any fixed k, 

Pk{p^u^) ^ Pkifiu) strongly in L°°(0, T; [F^ (0)]*) as iV ^ oo. (4.7) 

Now, we note that, for arbitrary k and N, and a.a. t G (0, T) there holds 



sup 



=1 Jn 



[p{p^u^{t))-Pk{p^u^)m^ 



sup 



p^'u^'m^-Pk^) 



=1 Jn 



< 



sup 



\ip - Pky^WL^ 



\m\tHi 



\ N N\\ 

\P U llL°°L2(t7) , 
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By a standard argument based on Rellich Lemma, one shows that 

sup \\lp - Pfc9j||L2(n) 

ll¥'llHi(n)=l 

as A; — )• oo. With the uniform bound on p^u^ in L°°(0, T; we can conclude that 

P^{p'^u'^) - P{p^u'^) strongly in L°°(0, T; as /c ^ +00, uniformly in N. 

(4.8) 

Of course, with a similar but simpler estimate, we also have 

P'^ipu) - P{pu) strongly in L°°(0, T; [//^(O)]*), as A: ^ +00. (4.9) 

Combining (14.71 ). (14.81 ) and (I4.9I ). we obtain finally that : P{p^u^) converges to P{pu) strongly 
in L^(0,r; [H^{Q)]*). Combining this strong convergence with the weak convergence of {u^) in 
L2(0, T; HI{Vl)), we might apply the method of P.L. Lions lEOl p.47] with the duahty bracket [HI{Vl)]*- 
H^{^) to prove that ^fp^u^ converges to ^fpu strongly in -L^((0, T) x Vt). Finally, we rewrite : 

I^{t)= [ y^u^ (S)^/^v^ -.Vej, 
Jn 

where : 

• \/p^u^ converges to strongly in L^((0, T) x Vt) 

• \fp^ converges to y/p strongly in L°°(0, T; L^i^)) 

• converges to v weakly in L'^{0, T; L^{il,)) (thanks to the imbedding H^{^}) C L^{Q)). 

Combining these statements, we get that In converges to / (with obvious notations) weakly in 
L\0,T). 

Such convergences result yield that {p'^, u^) satisfy c') for test functions tp of the form xi^)"^ 
with X £ ^([Oj^)) ^rid ijj G span({efc, k G N}). Via a classical density argument, the convergence 
extends to all 99 E i/^(0, T; n ^^(0, T; hI{VI)) such that 99|t=T = 0. 



4.4 Energy inequality 

We end this section by proving that the approximate solution (/o" , u") satisfies the further estimate : 



^ ' Jo Jn ^Pn Jo Jaf^ 

+ ^ /V - X z.|2 + n r / xSIn" - P^^n^p (4.10) 

^ps Jo Jds'^it) Jo Jn 

< I I p{-9)-u'' + [ pokol^ 
Jo Jn Jn 

for almost all t G [0, T]. For simplicity we drop exponent 7i in what follows. 

First, we note that the solutions {p^ , u^) of the Galerkin scheme satisfy (14.51) uniformly in N and 
that, up to the extraction of a subsequence \/p^u^ converges to in L^((0, T) x Vt). Hence, we 
may pass to the limit in (14.51 ) for almost all t G [0, T] . On the other hand, there holds: 
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By construction of the Galerkin scheme, Uq — )• uq in L'^{Q,) so that : 



n Jn 



lim / po\uq\'^ = I po\uo\'^ 

N 



Given the strong convergence of {u ) in L {{0, T) x $7) : 

lim r [ p{-g) .u'' = r [ p{-g) • n" 



Jn Jo Jn 



Given the weak convergence of in L'^{{0, T); (0)) and the strong convergence of in 
C([0, T]; LP(17)) we get that ^/Ji^D{u^) converges weakly to ^D{u) in ^^-^((O, T) x Q). 
In particular, in the weak limit, there holds : 



T r f-T 



< liminf /" I p^\D{u^^'^ 
Jo Jn 



10 Jn 

With similar arguments, we obtain : 

r [ Xs\u- Psu\' < liminf r [ - P^^^j^ . 

JO Jn Jo Jn 

Finally, we pass to the limit in the boundary terms. First, we introduce C/^ (resp. Ug) the 
extension Eq[u^] (the rigid vector field Pg u^) computed through the change of variable 0^q. 
As previously, we have: 

^ / |(n^-P5V)xH^ = r [ |(.^-nf )x.|2 = r [ \iU^-U^)x.\^ 
JdS'^it) Jo JdS^it) Jo JdSo 

Because of the weak convergence of and in L'^{0,T; H^^Q)), we have that and 
converge also weakly in L^(0, T; Hj-^^(M.^)) (see Lemma [17]) . Hence, we have also weak 
convergence of the traces on Sq. The lower semi-continuity of the L^-norm on dSo yields : 

^ [ \iu-Psu)xi^\^ = r [ \{U-Us)x,^\^ 

JdS{t) Jo JdSo 

< liminf r [ \iU^ - U^) X 
Jo JdSo 

< liminf/^/ \(u^ -P^u^)xiy\^ 



10 JdS^(t) 

Similar weak-convergence and semi-continuity arguments yield also : 

/ |n X i/p < liminf / / x i/p 

Jan Jo Jan 

This ends the proof of (I4.10I ). 
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5 Convergence 

In the previous section, we have obtained the existence of solutions of approximate fluid-soUd 
systems, namely satisfying a)-d). These solutions u"" are defined on some (uniform in n) time interval 
(0, T) such that 

dist(S""(t), > 25, for t G [0, T), for some fixed 5 > 0. (5.1) 

We must now study the asymptotics of as n goes to infinity, and recover a weak solution at the 
limit. 

In what follows, we will often make use of the notation 

{0\ := {x G M^ dist(2;,C') < r]} 

for O an open set and r] > 0. 

5.1 A priori bounds on n". Convergence in the transport equation 

The density clearly satisfies the uniform bound 

m.m{pF, ps) < < max{pF, ps)- (5.2) 
Combining (14.101 ) and (15.21 ) yields that 

ll^"lli«={o,T;L^(n)) + nWVTsiu"' - PSnn\\hmT).n) + II^^^K)lli^((o,T)xQ) < C. (5.3) 
for some constant C depending only on pp, ps, uq and T. 
In particular, up to a subsequence, the first inequality gives 

u" ^ n weakly* in L°°(0, T; Ll{n)). 

We can then pass to the limit of the transport equation d), using Proposition [5] The following conver- 
gence holds up to a subsequence: 

X§ ^ XS weakly * in L^{{0, T) x W"), strongly in C([0, T]; Lf„^(M3)) (p < oo), 

with 

Xs{t, ■) = ls{t), S{t) = (ptfiiSo) 
for an isometric propagator (f) = (pt^g G VF^'°°((0, T)^; C^^{R^)). Moreover, one has 
P^u"" ^ Psu weakly * in L°°(0,r;Q'^^(M3)), (^"^0 weakly * in Ty^'°^((0, r)^; ^^(M^)). 

In particular, one recovers the transport equation (12.31 ) setting us := Psu. 
Now, we can combine the second inequality in (15. 31 ). that yields 

Xsiu"" - Psu"") strongly in L^, 

with the strong (resp. weak) convergence of (resp. and PgU^). As n goes to infinity, we derive 
easily: 

Xs{u - us) = 0. (5.4) 
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Finally, the last inequality in (15.31) and Korn's inequality imply that 

We then introduce continuous extension operators 

En{t): {nGifi(F"(t)), divn = OinF"(t), n-i/|af^=0} ^ H^in), 

in the spirit of Corollary |7] As long as the S^{t) are 26 away from dQ, it is standard to construct these 
extension operators in such a way that 

\\En{t)\\cim)<Cs, VtG[0,r]. 

Hence, if we set Up{t, •) := En{t)u^ {t, •), we have that 

(m^) is bounded in L^{0, T; H^in)), (1 - Xs)iuF - u"") = 0, V n. 

From the ^^^(0, T; (0)) bound, we can assume up to another extraction that 

< ^ up weakly in ^^(0, T; Hl{n)). 

From above equality and from the strong convergence of x^, we then get: 

{l-Xs){uF-u) = 0. (5.5) 

Eventually, considering relations (15.41 ) and (15. 51 ). we get that the limit u of belongs to St- 
Hence, back to the definition of a weak solution, it only remains to show that the momentum equation 
(12.11 ) is satisfied by S{-),us,uf- 

5.2 A priori bounds on v". 

Prior to the analysis of the momentum equation, we must establish some refined estimates on the 
connecting velocity v"". We remind that was defined in Lagrangian like coordinates, see pai^agraph 
14.11 More precisely, 

v^t,^l,{y)) := dcl>l,\y(v'[U-it,.),U'Sit,-)]), 

where 

Enu^{t,^l,{y)) = d^l,\y{U^{t,y)), P^u^ {t,cPl,{y)) = dc^l.ly {U^{t,y)) , 

and = V^[U, Us] is some linear operator connecting U G H^{R^ \ Sq) to Us G H^CSq) over 
a band of width 6 outside Sq: see Corollary [8] We shall here specify our choice for the operator V^. 
Actually, we shall make it depend on n (V^ = y^'"), in order for the following additional assumption 
to be satisfied: 

\\V''"'[U,Us]-Uhr>aSo)s\So) < Cs,p {\m-Us)-u\\L,^QSo)+n'^'~'^^{U,Us)\\HH(So)s\So)) ' 

V2<p<6. (5.6) 



27 



We postpone the construction of such operator y^ " to the end of the paragraph. 
Back to f the additional assumption (15.61) impUes easily that for all 2 < ]5 < 6, 



(l-XS)K-n'^)||L^(0,T;L.(Q)) 

< Cs,p \\{u^ - PSumt, •) • + ^'/'"'/') • (5-7) 



But we know that 



(5"{t)) 



c 

< - (5.8) 

n 



where the last bound comes from the second inequality in (15.31 ). We emphasize here that, as the 
5"(t)'s ai^e all isometric to one another, the constant C does not depend on n, t. Interpolation with 
the similar other bound 



- PS'^"'){^1 ■) ■ 1 1 ^1/2 (g5n(t)) 



< CI ( ll'W"'||/fl((5n(t))^nJ;'n(t)) + ||u"||^2(5n(t)) ) 



(5.9) 



< C (ll'""'lljfi(F"(t)) + ll'""llL2(S"(t))) dt < C 
J 

yields that 

1^ \\{u^-PSu^){tr)-u\\l,^QS^^t))dt ^ 0, Vp<4. (5.10) 
Eventually, we get that 

||(l-xS)(t^"-^x")||i2(o,T;Lp(n))^0, Vp<4. (5.11) 
This will be much important in the treatment of the nonlinear terms. 



We conclude this paragraph with the construction of the operator V^'"' satisfying (15.61 ). We take 
U,Us in H^{R.^ \ So) X H^{So). Up to an extension of Us, there is no loss of generality assuming 
that Us G Hli 



Step 1. We shall construct a field V such that div V = 0, 

V\ds, = Us + {U- Us) -vv, and V\q(^s,), = U. (5.12) 

Therefore, we introduce a system of orthogonal curvilinear coordinates {si,S2,z) in a tubular neigh- 
borhood of OSq. Si, S2 are coordinates along the surface OSq, whereas z denotes a transverse coordi- 
nate. In particular, dS^ = {z = 0}. We set 

Id Id Id 

III osi 112 0S2 hz dz 
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the associated orthonormal vectors, with scale factors hi,h2,hz > 0. We remind that 

V/ = ^dsje, + ^dsje2 + ^dje, (5.13) 
hi /i2 hz 

for a scalar function /, whereas 

/ = T^T- idsAh2hJi) + ds,{hhj2) + d^ihM,)) (5.14) 
nin2iiz 

for any field f = fiei + f2e2 + fzSz- We then set 

Vi := (1 - x{nz)) U + x{nz){Us + [{U - Us) ■ e,] e,) 



for a smooth truncation function x : — ^ [0, 1] equal to 1 in a neighborhood of 0. Clearly, for all 
p<6,andl + i = i, 

ll^i-f/||L.((5o)A5o) < C,,sn^'/^'\\{U,Us)hHiSo)s\So) (5^5^ 

< C'p,5n^'^'^m,Us)\\HmSo)s\So)- 
Also, Vi satisfies (15.121 ). But it is not divergence-free: formula (15.141 ) yields 

C 

div Vi = x(nz)div {[{U-Us)-ez]ez), so that ||div 14||l2((So)A5o) ^ -^\P-^s\\maso)s\So)- 

y n 

To obtain a divergence-free field, we note that both U and Us have zero flux through dSo and use 
Proposition [6l there exists a field V2 such that 

div V2 = -div Vi in (5o)5 \ So, V2\aso = ^2\d(So)s = 0' 

and 

II^2||hi((So)ASo) ^ ^ 11^- f^5||Hi((5o)A5o)- 



In particular, by Sobolev imbedding, one has for all p < 6 



II^2||lp((5o)A5o) ^ -7^\\(U^Us)\\m{{So)s\So)- (5-16) 



Finally, the field V := Vi + V2 fulfills our requirements. 
Step 2. We construct a field W such that div W = 0, 

W\aSo = m - Us) • ly] u, and l^|a(So), = 0- (5.17) 
In the same spirit as in the first step, we take 

Wi := x(y)[(C/-[/s)-H.=o]e. 

where x is again a truncation function: x = 1 0' ^rid x = outside [—1,1]. A rapid computation 
shows that 

\\Wi\\LHiSo),\So) < CsWiU - Us) ■ i^hpiaso), Vp (5.18) 
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By Proposition |6j there exists a field W2 such that 

div W2 = -div Wi in {80)5 \ Sq, W2\dSo = ^'2\d{So)s = 0' 

and 

ll^2||//i({So)A5o) ^ \\{U - Us) ■ i^hHaso)- 
In particular, by Sobolev imbedding, one has for all p < 6 

\\W2\\Ln{So)s\So) < Cs \\{U - Us) • lyWmaso)- (5-19) 
Finally, the field W := Wi + W2 fulfills our requirements. 

Eventually, we set 

'V^[U,Us] = [/ outside {So)5, 
< V^[U,Us] = V-Win {So)s\So, 
.V^[U,Us] = Us in So. 
Combining KTSl . KTB . (I5J8]) and (l5J9l ) leads to K6\) . 

5.3 Approximation of the test functions 

The weak formulation of the momentum equation involves discontinuous test functions 99 G Tr^ 

^={l-Xs)ipF + XsVS, eV{[0,T);V^m), V95GP([0,T);7^), 

with 

On the contrary, the approximate momentum equation c) involves continuous (or at least H^) test 
functions. Hence, we will have to approach 99 by a sequence (99") in L^(0, T; H^{Q,)). Due to the 
discontinuity of the limit, the •)'s will exhibit strong gradients near dS"'{t). Precise estimates 
are needed, that are the purpose of 

Proposition 12 Let a > 0. There exists a sequence in VFi'°°(0, T; (f))) nL°°(0, T; HI{Vl)), 
of the form 

99" = (1 - + xgv'S, 

that satisfies 

• WVTsi^'h - Vs)\\c{%nLnn)) = 0{n~^/P)forallp G [2,6]. 
, (^ri ^ ^ Strongly in C([0, T]-L^{Q.)). 

• \W\\c{[o,T];mm = 0{n'^/'^). 

• Wx'hidt + • V) (99" - (^5) llL=°(o,T;L6(f7)) = 0(n-"/6). 

• (5t + P^u" • V)(^'^ ^ (Qt + Psu ■ V)(^ wea% * in L°°{0, T; L6(J7)). 
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Proof of the proposition. 

The point is to build a good approximation 99^ of 995 over the solid domain. Broadly, we want 

V'sl*' ■)las'"(t) = '^F{t,-)\Qsr^{t) Vi, 

and 

•) w V95(t, •) in away from a n"" neighborhood of ^^"(t) Vt. 

Therefore, we proceed as for v", by using lagrangian coordinates: we define and through the 
formulas 

{t, <Ploiy)) = d<plo\y (<f y)) , {t, ^loiy)) = dcp^^ly {^pit, y)) , 
and the goal is to define properly some related to ip'g by the formula 

Note that $5 and depend on n through the propagator <j)^, but we omit it from our notations. The 
only thing we have to keep in mind is that the bounds on guai^antee that and <I>i? are uniformly 
bounded in W^'°°{0, T; Hj'^^iR^)) for all k. 

Thanks to the change of coordinates, the problem is now in the fixed domain Sq. Roughly, we 
want to build in such a way that 

nit,-)\dSo='^F{t,-)\aso Vt, 

and 

^s{t, •) w •) in 5o away from a n~" neighborhood of dSo V t. 

Note that time is only a parameter in the system. The construction of follows the one of V, 
performed in the previous paragraph. Step 1. We take <I>g under the form 

The first term has the explicit form 

^5,1 = + ii^F - ^s) - [(^s - ^f) • e.] e,) . 

Again, x is a smooth truncation function near 0, and 2; is a coordinate transverse to the boundary: 
dSo = {z = 0}. It is easily seen that <&g ^ satisfies the right boundary condition at OSq. Moreover, 

11^5,1 -^s||iyi>-(0,T:Lf{So)) < Cn-'^/P Vp<oo, \\^s,i - ^s\\w^.^{Q,T;m(So)) < C n""/^ . 

(5.20) 

But it is not divergence-free. By applying formula (I5.14I ). we get 

div = x(n° z)f, f := div (x ^ {{<!>p - <^s) - [{^S - • e,] e,)) . 

In particular, j" is uniformly bounded in W^'°°{0, T; L^(5o)) 
By Proposition |6l there exists some field $g 2 satisfying 

div ^>?2 = -div «>gi in ^o, '^s,2\dSo = 0, 
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and 

11*^5,2 lliyi'°°(0,T;Hi(5o)) - llx("-"^)j"llH/i.°°(0,r;L2(5o)) < Cn"""/^. (5.21) 

In particular, 

11*^5,2 II w^i'°°(o,T;L6(So)) - Cn~'^/'^. (5.22) 

Back to the moving domain (in variable x), one can combine the estimates (I5.20I) - (I5.21I) - (I5.22I) 
with the uniform bound on 0" in Ty^'°°(0, T; C°°{Q)). From there, one can deduce the estimates of 
the proposition. For the sake of brevity, we only treat the two last items. Namely, we write 

d 

WX^sidt + PsV^ ■ V) (<^" - ^s) llL-((0,r);L6(n)) < C ||^ ^^^oly (^S " IIl-((0,T);L6(So)) 

where the last inequality involves (15.201 ) and (15.221 ). This bound implies in turn that 

{dt + P^u^-V)^^ = (l-x§)(5t + P>"-V)(^i. + x§(5i + P>"-V)(^5 + 0(n-"/6) m L\n) 

The products at the r.h.s. are then easily handled using the strong convergence of x§ (and the weak 
convergence of P^n"). We obtain 

{dt + P^n" • V)99" ^ {dt + Psu ■ V)(/? weakly *in L°°(0, T; L^(Q)) 

as expected. This concludes the proof of the proposition. 

5.4 Convergence in the momentum equation: linear terms 

We now have all the elements to study the asymptotics of the approximate momentum equation c). 
Given an arbitrary 93 G 7t> we consider an approximate sequence (99") as in Proposition [T2l We shall 
take as a test function in c), and let n tend to infinity, so as to recover (12.11 ). We shall rely on the 
fields Up and u f introduced in paragraph 15.11 We remind that 

{1-Xs)uf = (I-Xg)^", ^ UF weakly in L\0,T; Hl{n)). (5.23) 

To lighten notations, we shall write := P^u^, us := Psu. We remind that these rigid fields 
satisfy 

u'^s ^ us weakly * in L°°(0,r;iy;^^(]R3)) VA:. (5.24) 

In this paragraph, we consider the asymptotics of all terms but the convection one. 
• We write the diffusion term as: 



Jo Jn ^ Jo Jn 



From the strong convergence of to in C([0, T]; LP(0)), and the weak convergence of 
to UF in L'^{0,T;H'^{n)), we deduce 



/r ^ /" f 2iiF{l-Xs)D{uF):D{ipF). 
Jo Jn 
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As regards > we use the bounds 

\\^Diu^)\\hmT)xn) = 0(1), y^L^io,T;HHn)) = 0{n^'^) 
established in the previous pai^agraphs. They imply 

If we choose a <2, then I2 goes to zero as n goes to infinity, and finally 

r [ 2/."Z)(n") : Di^^) ^ T [ 2/i,.Z)(nF) : Diipp). 
Jo Jn Jo Jfu) 



The boundary term at dO. reads 



1 I f 1 

(n" X v) ■ {ip"' X ly) = — — / / {up x v) ■ {ipp x i^) 



2/3c Jo Jdn 2/?Q 

1 



2/3q 

by the weak convergence of in L^(0, T; {^))- 



Jdn 

T f 

/ {up X v) ■ i^F X l^) 

Jdn 



We deal with the boundary term at dS"" as in the Galerkin approximation. We introduce := 
PgU^, := Pgip"^ = Ps'f's ^'^'^ capital letters to denote velocity fields when seen through 
the change of variable. We then have, as in the computation for the Galerkin method : 



1 



((n^ - ng) X 1.) • ((v.^ - rg) X z.) 



2/3s Jo JdS"{t) 
1 



"2/55 



[ ((C/]^-C/^)x^).(($p-i?g)x^), 

JdSo 



where we used that c/?^ = (pp- We note here that 99" converges to 99 in C([0, T]; L^(il)) so 
that combining with the strong convergence of Xs it yields that rg converges to rg := P^t^ in 
L2(0, T; HI^^{R^)). Through the change of variable. Lemma [17] yields that : 

i?g ^ Rs strongly in L^{0, T; H^/'^{dSo)). 

Then, we combine the respective convergences of u"', Ug with LemmafTTlvielding. with obvious 
notations : 

Up weakly in L'^{0, T; H\n)), Us weakly in ^^(0, T; H^{n)) . 

We apply these convergences together with the continuity of traces on dSo CC 0,, and go back 
to the moving geometry, to obtain finally : 



1 



2/?S Jo JdS"{t) 



((u" - ng) X u) . ((^" - (^g) X ^) 



T 



2/^5 Jo JdS(t) 



{{up - Us) X u) ■ {{ipp - LPs) X u) 
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to treat the penalization term we use the bounds 

^llv^(^"-^>")lli2((o,T)xQ) =0(1), \\V^si^'s-^s)\\c{[m-,L^m=0{n-^l^) 

estabUshed in the previous paragraph. We also remind that ips, as a rigid vector field, satisfies 
(ps = Ps^s- From there, 



Jn 



Jn 
Jn 



If we choose a > I (which is compatible with the former constraint a < 2), the penalization 
term vanishes as n — )• +oo. 



5.5 Strong convergence of 

To show that {S, u) is a weak solution over (0, T), we still have to pass to the limit in the convection 
term ^ 

conv" := / / /)" (n'^ • dt^'' + i;" (g) : V(/?") . 
Jo Jn 

To compute this limit, we first prove 

Proposition 13 Up to the extraction of a subsequence, (u") converges to u in L^((0, T) x 0). 

This result is obtained applying the method introduced in the reference Il23l (see also ifTOl for the 3D 
case). We first introduce some notations. Given < s < 1 and S a bounded connected subset d f], 
we denote 

n'[S] = the closure of {v G hI{9.) such that G U] mH'{VL). 

As TZ'^IS] is a closed subspace of H^{Q,) we denote P^[S] the orthogonal projector from H^{Q) 
onto this subspace. Given s' > s, we recall that TZ^ [S] is a dense subspace of 7^''[S'] , and that the 
imbbeding TZ'^ [S] C Tl^[S] is compact. If s = 0, we shall drop exponent s. We emphasize that in the 
case s = the projector P[S] does not coincide with the Ps introduced in (13.11 ). 

Our first step is the following approximation lemma : 
Lemma 14 Let s < |- 

i) The sequence {un) is uniformly bounded in L^(0, T; H^(Q,)). Moreover, there is e = e{s) > 
such that for all h < 6/2, 







(5.25) 



ii) One has u G L^{0, T; H^[Q)). Moreover, there exists e = e{s) such that for all h < 6/2, 







\u{t,.)-P'[{Sm]u{trWHs^^-^dt < Ch' 



(5.26) 
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where, in both cases, the constant C depends only on initial data. 

Proof of the lemma. We only prove the first item of the lemma, the second one being simpler. It 
relies on the construction of a suitable approximation v"^ of n", rigid in a /i-neighborhood of 5". This 
approximation will satisfy the following properties : 

• {v"^) is bounded in L^(0, T; H^{Vl)) for s small enough 

• vl{t, ■) = Psnu'^it, •) in and v'/iit, •) = n"(t, •) outside iS{t))s for a.a. t G (0, T). 
Note that it implies v]^{t, ■) G 7^*[(5'"(^))/,] for a.a. t £ (0, T). 

• for h sufficiently small and for a.a. t G (0, T) there holds 

\\u^{t, •) - vUt, ■)\\LHn\iSHt)h) <Ch\ {mn^it. OIIl^cq) + OIIh^f-W)) 

+ C\\{u^-P^u^)-vU2^esHt))^ 
II vl{t, Olk^co < C(l + h^) (||P^n"(t, OIIl^c^^) + lk"(t, OkMi^-W)^ 
+ C||K-P>-)-i.|U2(95„(i)). 

(5.27) 

Before giving further details on the construction of v"^ we explain how the previous properties imply 
Lemma l5.25l By interpolation of (15.81) and ( 15.91 ). we obtain 

UK - Psy^na. •) • HhidS-it)) dt < (5.28) 



n 

We square the inequalities in (15.271) and integrate from to T. Using (15.281) with the uniform bounds 
(15.31 ). we end up with 

rT \ 1/2 

(5.29) 



Moreover, 



Ih" - ^^hlli2((5"(t))h) ) =( / 11^" --fs^"llL2{(5"{t));,) 



< ( / Ik" - -Ps^"llL2(S"(i)) ) + ( / II'"" ~ -fs'^'^llz,2((5n(i))^^\5n(i)) 



T \ 1/2 / /.T X 1/2 



Using (15.31) . we get 

..n pn n\\2 

^ ^S^ \\L-i{{S^{t))h\S'^{t)) 



< ^ + cVh( [ 

< + CVh ( f llu" - P?n" 



1/2 



/ .T xl/2 ^ 

vn 
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Combining this last inequality with the first inequality in ( I5.29I ) yields 

Ik" - ^^hllL2((o,T)xn) < C {h"'' + (5.30) 
As regards the norm, s < 2/5, we use the second inequality in ( I5.29I ) to write 

- K\\L2{0,T;H^{n)) < ij^ Ik" - -Ps^"'llff''(S"(t))^0 ~^ ( /g ll'""llH''(F"(t))'^M 

/ .T X 1/2 

+ n \K\\Hs(FHt))dt] < \K - Psul\L^iO,T;H'm +0(1) 



Finally, we have 

Ik --Ts^i llL2(0,T;H-(n)) < 11^^ - ^S^ llL2{(0,T)xt7) Ik "^5^ IIl2(0,T;//1(!^)) 

, 1 1 

< C —p= < C as soon as s < 



n J 3 



We end up with 

IW^ -Vh\\L^{o,T;H''{n)) < C as soon as s<-. (5.31) 
One last interpolation between (15.301 ) and (15.311 ) shows that for all s < 1/3 and e = e(s) > 0, 



\u 



- Vh\\L^iO,T;HHn)) < C{h'+n 



As vj^{t, •) belongs to 7^**[(5"(t))/j] for all t, by definition of the projection, the same inequality holds 
replacing v]^ by P[{S"'{t))h], as expected. 



We still have to achieve the construction of v^. It follows the construction of v"', cf paragraph 
15.21 It is actually simpler, because we only look for a vj^ with H'^ regularity for small s. In particular, 
jump on the tangential part at d{S")h and d{S^)s will be allowed. 

As before, we go back to Lagrangian coordinates : we look for a vj^ under the form 

v]i{t,(i>tAy)) = d4>t,o\yV;:{t,y). 

Also, we define C/" and through 

In this way, we are back to a static problem. For brevity, we shall omit temporarily the time dependence 
in our notations. The point is to build a field Vf^ satisfying 

= in {So)h: Vf: = [/" outside (5o)5, 

and suitable estimates. 

Therefore, we follow paragraph l5.2[ We parametrize (S'o)5\5o by curvilinear coordinates {si,S2,z), 
z being the distance at dS^. Hence, d{So)h = {z = h}. Then, we introduce 

:= {l-X (^)) + X (^) (U^ + Kf/" - U^) ■ e.] e.) 
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and the solution of 

r = -divF^I^i, miSo)s\{So)h, 

\ V^^s = on 5(5o)5 and 9(5o)h 

Computations similar to those of pai'agraph I5.2l yield : 

IIK^i - t^"llL2({So)A{So)h) ^ h-i\\{U'',U^)\\Hi((So)6\{So)h)x'R-, (5-32) 
IIK!lll//i({So)A(So)h) ^ ^~^||(^^",f/5)ll/fi((5o)A(5o)h)x7^> (5-33) 

\\Vh,2\\Hmso)6\{So)h) ^ Ch2\\{U'^,U^)\\Hi((^So)s\{So)H)xn- (5-34) 

Let us emphasize that the constant C in the last inequality can be chosen uniformly in h, see lITTl 
Theorem III.3.1]. It follows by interpolation that 

2 — 5s 

\\K,i + K,2\\H-aso)s\{So)h) ^ C'^" ll(f^"'f^5)ll^^l((5o)A(5o)h)x7^• (5-35) 
Finally, we build some WJ^ = VY"^ where Y"^ is the unique solution of : 

Ay^" = in (5o)A (5o)h , . 

^zY;: = (C/^ -[/")• , ond{So)h, such that / 1^ = 0. 

a^y" = 0, ond{So)5, JiSo)s\{So)h 

we recall that = on d{So)h- By standard elliptic regularity results, there exists a constant C 
independent of h such that : 

\\Wh\\LmSo)s\{So)H) ^ \\Yh\\HmSo)s\(So)h) ^ - C/") • e^||^-i/2(a(5(,)^) , 

By interpolation, we get 

II^."IIhV2((5o)A(5o).) < - Un ■ e.WLHdis,),) (5.36) 

Now, we write 

\m-un-eA\L^i9i^s,),) 

< C {hHviU^s - f^'^)llL^((So)A(So).) + \m - Un ■ 6.11^2(950)) • (5-37) 

Eventually, we set Vj^ := V^^i + Vi^2 ~ ^h- We stress that the normal component of Vfi" is 
continuous across d{So)h and d{So)s- Hence, for any s < |, the H'^ norm of Vf^ over the whole 
domain is controlled by the sum of the norms over {So)h, {So)s \ {So)h and Q"" \ (50)5, where 
Q"^ is a shorthand for cI)q^{^). It follows from this remai'k and the previous inequalities that: for all 

s < 1/2 

Also, one has 

- ^"llL2{C"\(5o)h) 

< C (h'/'{\\U^H^^p.) + \\U^\\n) + \\m-Un-e.h2i9So))- (539) 
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Back to the moving domain, and accouting for time dependence, we obtain (15.271) . 



Tlie second step in tlie treatment of tiie nonlinear- terms is a control of the Hausdorff distance 
between S"'{t) and S{t') for close times t, t' G [0, T]. This is the purpose of 

Lemma 15 Let h > 0. 

i) There exists no > such that for all n > tiq, 

S^{t) C (S(t)),/4 C (5"(t))v2 yt€[0,T]. 

ii) There exists ij > such that for all tQ £ [0, T], for all t £ [to — r], to + ijjn [0, T] 

{S{t)\/^ c {S{t,))n c {S{t))2h. 

Note that condition (15.11 ) and point i) of the lemma (applied with h = 6) imply that 

3 

dist{S{t),dn) > -5, for t e [0,T], for some fixed 6 > 0. (5.40) 

Proof of the lemma. We first treat i), focusing on the first inclusion (the other one is proved in 
the same way). To this end, we recall that the associated sequence of characteristic functions 
converges to X5 in C([0, T]; (0)). This implies that 

sup |5"(t) A5(t)| = sup ||x5(i,-) -X5(t,-)llLi(f7) ^Owhenn^oo, 
te[o,r] te[o,T] 

where we denoted A the symmetric difference of subsets of M'^. Let us now take /i > and assume a 
contrario that there exists a sequence of time tjt G [0, T] and of integers nj^ going to infinity such that 

As 5(tfc) is isometric to 5o, which satisfies: 

3 r > s.t. for all x £ Sq there exists a euclidean ball B with radius r satisfying x £ B d Sq 
there exists for all k a ball i?^ with radius r' = min(r, /i/16) such that 

B'^ClS^-{tu)\S{tu), 

so that 

sup A5(t)| >^^, 

tG[0,T] <3 

which yields a contradiction. Consequently, there exists no such that, for all n > no, 

s^{t) c (5(t))^/4, vtG [o,r]. 

The second item ii) is obtained in the same way. Let /i > and assume for instance that the first 
inclusion does not hold. Arguing as previously, we are able to construct two sequences {t^) and (t^) 
converging both to G [0, T] and such that S{t^) \ S{t^) contains a ball of fixed radius. Once again, 
this contradicts the continuity in of xs at to. 
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Thanks to the previous lemmas, we can conclude the proof of Proposition [T3j following very 
closely Il23l . At first, very minor adaptation of the proof of 1231 Proposition 7.1] yields: for s G 
(0, 1 /3), there exists ho, such that, for all h € (0, Hq): 

lim r [ p'^u'^ ■ P'[iSit))hK = r [ pu-P-'[iSit))h]n. (5.41) 
"^°<=yo Jn Jo Jn 

We remind that the main idea behind this limit is the following: thanks to Lemma [T5l for any field 
^, the projected field P'^[{S{t))h]{S,) is rigid in a neighborhood of S"'{t) for n large enough. Hence, 
if one uses P'*[(S'(t))/i](^) as a test function in the momentum equation, the boundary term at dS^{t) 
and the penalization term vanish: roughly, one recovers a uniform bound on (9tP**[(5(t))/i](/>"n") in 
a Sobolev space of negative index, and from there compactness. For all details, see ll23l Proposition 
7.1]. 

Then, one establishes that 



lim 



Jn 



Jn 



(5.42) 



^P" = - Xs) + PsXs^ P = - Xs) + PSXS^ ■ The idea is to write 

V P"KI^- r [ P\n\' 
Jn Jo Jn 

^ / pV-P^[{Sm]iu-)- r [ pu.P^[iSit)),]iu) 
Jn Jo Jn J 

+ r [ p^'u^ ■ - p'[s{t)hK)dt + r [ pu-{p'[s{t)h]u-u)dt 

Jo Jn Jo Jn 



The first term at the r.h.s. is controlled using (I5.41I ). whereas the last two are treated thanks to Lemma 
[T3l note that {S{t))h C (5'"(i))2h for n large enough by Lemma [T5l so that 



JO 



The final step of the proof consists in showing that 

rT r r-T 



Jn 



p\u 



n\2 



Jn 



p\u\ 



which yields the strong compactness of uP" {p having positive lower and upper bounds). The idea here 
is to write 



T 



Jn 



p{\u 



n,2 |2. 



< 



Jn 



{pVI'-pH') 



+ 



Jn 



ip--p)\u 



n\2 



The first term at the r.h.s. goes to zero by (15.421 ). For the second one, we use that p" — )• p strongly in 
C([9, T]; L'P{Q)) for all finite p and that is uniformly bounded in for some p' > 1, thanks to 
the uniform bound on n". Again, we refer to ll23l for all details. 
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5.6 Convergence in the momentum equation: nonlinear terms 

Thanks to the strong convergence of Proposition [131 we are now able to split conv" in a suitable way. 
Let us first remind that v"" is identically equal to m§ inside S^, whereas is identically equal to ifp 
outside S*". This allows us to decompose the convection term as follows: 

conv-= r [ pFil-Xs)uF-dtV^F 
Jo Jn 

+ [ [ Pf{1 - Xs^"" ®u],:V^F+ [ psxm + • V)(^" • n" := If + + 1^ 
Jo Jn Jo 

The convergence of /" is cleai" 

^ [ [ Pf{1 - Xs) up ■ dt^PF- 
Jo Jn 

The convergence of follows from the fourth item in Proposition [T2l which cleai^ly implies that 



Is 



3 = r [ PsXm + ■ V)y^5 • + 0(1). 
Jo Jn 



Using the strong convergence of x^""" to xs^s in ^ ((0, T) x il), it is then easily shown that 

1^=1 [ PsXsOt^Ps ■ us + [ [ PsXs^s-^'Ps-u'^ + o(l). 
Jo Jn Jo Jn 

Now, we write the second term at the r.h.s. as 

[ [ psXsUs ■ V</9s -u"^ = [ [ PsXsUs : Vv?5 
Jo Jn Jo Jn 

I psxsus ® «s : Di^s) = 
Jn 

as ips is a rigid vector field. 

It remains to study • We know from paragraph 15 . 21 that 

(1 - - n") = (1 - - u'^) ^ Oin L\0,T; LP{n)), Vp < 6. 

It follows that 

IS = r [ Pf{1 - Xs)uf ® «F : Vyp^ + o(l). 
Jo Jn 

= [ /pi.(l-X5K^«": V(^i. + o(l). 
Jo Jn 

In this last identity we collect the strong convergences of to u in L^((0, T) x il) and of Xg to x in 
C([0, T];L^^{n)), together with the uniform regularity of (n", ti) in L^{0, T; H^/^{Vl)) (see Lemma 
[141 ). which yields that {un,u) are uniformly bounded in L^(0, T; L^'^/^^(0)). We obtain then : 



lim / / pf{1-Xs)u'' ^u"" -y^F = f f pF{l-Xs)u^u:VipF 

>-^°°Jo Jn Jo Jn 

= pp{l - xs)uF <8) up : V(pF 

Jo Jn 



This concludes our proof. 
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5.7 Energy inequality and extension to collision time. 

We pass to the weak limit in (14.101 ) and prove that the solution (p, u) satisfies the further energy 
estimate (I2.2I ). First, we note that (14.10b implies 

Jo Jn ^Pn Jo Jan 



+ ^ r / - ^^«") X ^1' < r / Pi-9) ■u'^ + [ Po \uo 

^Ps Jo Jds'^it) Jo Jn Jn 



2 

MO I "0 1 

n 



for all n. As we have convergence of ^J^u^ in L^~^((0, T) x $7) we can pass to the weak limit in 
this inequality for almost all t G [0, T]. As S(i) remains far from 517 we treat boundary terms in a 
similar way as in paragraph 14.41 The only term which requires a new treatment is : 

2;u"|Z)(u")|2. 



Jn 



For this term, we note that, because of Lemma [TSl there holds for arbitrary /i > and n sufficiently 
large : 

Jo Jn\{s{t))h Jo Jn 

If we let n go to infinity, and then h go to 0, we obtain : 



/ / 2/iir|L»(u^)p < liminf /" f 2fi'''\D[ 
Jo Jpit) Jo Jn 



u 



lF{t) 

for almost all t G [0, T] . Hence, passing to the limit in (14.101) yields (12.21 ). 

Our solutions are limited in time to avoid collision. Namely, the only shortcoming of our con- 
struction is that it requires the distance between S{t) and dil to be larger than a fixed positive distance 
6 through time. However, as long as we are given an initial data uq G L'^{Q) and an initial position 5*0 
such that 5o (s we are able to construct a small time T depending only on the inital position of 
in and the norm of uq such that the solution exists and satisfies this property on [0, T]. As our 
solutions satisfy also energy estimate (12.21 ) we might reproduce the arguments of lITOl Lemma 2.2] to 
concatenate solutions in time and prove existence of at least one weak solution until collision time. 
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A Weak/Strong convergence and isometrics 

In this appendix, we study the influence of isometric transformations on weak and strong convergence 
of sequences. First, we prove : 

Lemma 16 Let E C([0, T];Isom{R^)). Given w : (0, T) x ^ j^^^^ 

w{t,cl){t,y)) ■.= dcl)t\yWit,y), V (t, y) G (0, T) x M^. (A.l) 

Tlwn 

• Ifw G L2(0,T;F1(M3)) then W £ L'^{0,T; (R^)). 

• Ifw G C{[0,T];H^{R^)) then W G C{[0,T]; H^{R^)). 

• The same assertions hold true replacing H^{R^) by Hl^^(R^^). 

The proof of this lemma is based on the fact that formula (lA.ll) for fixed t defines a unitary transfor- 
mation of H^{R^). The details are left to the reader. Second, we obtain : 

Lemma 17 Let 4>^ : [0,r] x such that (j)^{t, •) G Isom{R^) for all t G [0,r]. We assume that 
(p^ converges to (j) in C([0, T]; C^^{R:^)) . Given a sequence (w^) : (0, T) x M'^ — )• R^ we define : 

u;^(t,<A^(t,2/)) :=#f|,W^^(t,y). 

Then, with obvious notations: 
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• If (w^) converges to w strongly (resp. weakly) in L^{0,T; H^{M.^)) then (W^) converges to 
W strongly (resp. weakly) in L'^{0, T; H^{M.^)). 

• If{w'^) converges to w in C([0, T];H^{M.^)) then (W^) converges to W in C([0, T]; //^(M^)). 

• The same assertions hold true replacing i7^(M'^) by Hl^^(K^). 

Remark. We point out that and satisfy symmetric relations: 

so that fields ]¥"■ and w^, resp. W and w can be switched in this lemma. 

Proof of Lemma\17\ We first remind that (j)^ is an affine isometry, so that (for all N, t) 

\d(t)f\yx\ = \x\ , \[d(t>f \y\-^ M d(t>f \y\ = |M|, y{x,y) G x VM G M3(R). (A.2) 

The same relations hold for cj) instead of (p^ . 

Strong convergence. We focus on convergence in C([0, T]; i/^(]R^)), the strong convergence in 
L^H^ being ti'eated in the same way. First, we note that our previous lemma yields: 

G C([0,T];ifi(M3)) for any Af, G C([0, T]; /^^(M^)). 

Then, we write 

II^^^-^^IIC([0,T];L2(M3)) < sup If (t) + I2 {t) + SUp /f (t) 

where 

|lf(t)|2:= / \w''{t,y)-d^f o[d^^]-\W{tMt\-^°<i^f{y))\^dy, 
\l2{t)? := / I [#f o d4>t\y W{t, [4>t]-^ o 4>f{y)) - W{t, [cPt]-' o 0f (y))|' dy, 

and 

l^"" / \W{tM-'o^f{y))-W{t,y)fdy. 
Using (IA.2I ). we have easily 

\Iiit)\' = [ \w''{t,<pf'{y))-w{<P^{y))fdydt 

[t, x) — w{t, x)]^ dx, 

which tends uniformly to when — )• oo by assumption. We then get 

\I^it)\^ < SUp|d<iOd,/.t,o-/d|' / \Wit,[^t]~'o^f{y)fdy 
t,y JRS 

< sup \d(t)Q^t o dcptfl - Id\^ ||W^||c([0,T];//i(R3)) 0. 

t,y 
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Finally, the continuity of W with values in L^(M^) implies that: 

\W{t,y)\^dy 



L 



\y\>A 

can be made arbitrary small uniformly in time, taking A sufficiently large. So, we apply the local 
convergence of (p^ to 4> to obtain that, for N sufficiently large, there holds : 

/ V 1/2 / X 1/2 

\Iz{t)\ < I \W{tMtr'o4>^{y))\^dy] + / \W{t,y)\^dy] 
\J\y\>A / \J\y\>A / 



+ ( / \W{tMt]''o'P?{y))-W{t,y)\\y\ 
\y\<A / 







)■'■* ( 


7 

J\y\<A 



<2[ I I |H^(t,(y))n + ( / \W{t,[<Ptr'oct>^'{y))-W{t,y)fdy 



■J\y\>A/2 

The first term at the r.h.s. is independent of N and goes to zero as A goes to infinity. Moreover, for 
fixed A, [(t>t\~^ ° 4>t^{y) converges to y uniformly in [0,r] x {\y\ < A}. Hence, for fixed A, the 
second term at the r.h.s. converges to zero as N goes to infinity (continuity of translations in L^) 
uniformly in t. We conclude that goes to zero, so that converges to in C([0, T]; L^' 
The convergence of VW^ to VW follows the same lines, which yields the result. 

Weak convergence. Again, we only prove the convergence on . The convergence in Hj^^ 
is similar. We assume here that {w^) converges to w in L^(0, T; iJ^(M'^)) weak. Given x ^ 
C;?°((0, T) X M3) there holds : 

^ f W^it,y)-xit,y)dtdy = T [ w"" {t,4>? {y)) ■ dct>^,^\yx{t,y) dt dy 

..N u \ /jr^A^i-li ..u r^A^i-l 



w''{t,x)-[d[<f>n-\Yx{tA<Ptr{x))dtdy, 

where we applied again that d4)^\y is a hnear isometry. Because of the strong convergence of in 
C7([0, T]; Qoc(M3)) there holds : 

{dm-%y\{t, [./-n-^lx))^ (dM-iU.)"'x(t, strongly in L2((0,r) xM3) 

so that, combining with the weak convergence of we obtain : 

W''it,y)-x{t,y)dtdy^ r [ W{t,y) ■ x{t,y) dt dy . 





Similar arguments yield also that: 



[ [ VW''{t,y):Vx{t,y)dtdy^[ [ VW{t,y) : Vx{t,y) dt dy . 

Jo Jo 

which ends the proof. 
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